Molecular dynamics simulation of polymer insertion into lipid bilayers 
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Relatively short peptides, such as toxins and antimicrobial-peptides, are known to insert them- 
selves into cell membranes. On the basis of simple bead-spring models for the membrane lipids, the 
peptide, and water, detailed processes of the peptide insertion is investigated by molecular dynamics 
simulation; our special concern is in the highly cooperative motions of membrane lipids and the pep- 
tide. Our model lipid has a head group of three hydrophilic beads and a tail of seven hydrophobic 
beads, while the model peptide is a block-copolymer made of hydrophilic and hydrophobic blocks, 
with total length of 200 beads. In addition, each water molecule is represented by a single bead 
which has considerably larger interaction energy. We first confirm that our present lipid model can 
support spontaneous formation of bilayers in water. Then we place the model peptide near the 
bilayers and monitor the microscopic process of the adsorption, insertion, and translocation of the 
peptide. When the peptide molecule is set free to interact with the membrane, the hydrophobic 
blocks of the peptide strongly favor intimate contact with the lipid tails. This strong attractive 
interaction gives rise to a severe membrane perturbation, which leads to the formation of a pore, 
though short-lived, in the membrane. It is found that the side-surface of the pore is almost covered 
with the hydrophilic heads of the lipids, which seems to help the hydrophilic blocks of the peptide 
translocate across the membrane. We also monitor the rate of flip-flop inversion of the directions 
of lipids, and we find it markedly increases during the peptide contact and insertion. The peptide 
insertion/translocation, the formation of membrane pore, and the flip-flop motions of the lipids are 
thus found to be closely interconnected. 

PACS numbers: 



I. 



INTRODUCTION 



Lipid bilayers are common frameworks of various bi- 
ological membranes. Microscopic structures of the bi- 
layers, molecular level understanding of their formation, 
various morphologies they exhibit, and their interactions 
or fusions, are all of fundamental importance in the stud- 
ies of cell membranes [1]. Most of the exquisite functions 
of the cell membranes are, however, carried out by coop- 
erative work of lipids and membrane proteins: receptors, 
ion channels, etc. In fact, proteins are actively moving 
even within the cell by inserting into or translocating 
across membrane systems between various organelles [2] . 

Peptides and proteins usually accomplish insertion or 
translocation either assisted by complex proteinacious 
machinery or through fusion of membranes [2]. Rela- 
tively small peptides such as toxins and antimicrobial 
peptides, however, are known to insert themselves by 
physicochemical processes, in which four- stage or two- 
stage models have been proposed for the insertion and 
reorganization [3-5]. The peptides first get adsorbed to 
the membrane surface where they change their conforma- 
tion into a-helix due to large hydrophobic environment 
at the surface. Then they start insertion, whose molecu- 
lar mechanism is least understood, followed by aggrega- 
tion into an organized structure. The driving force easily 
think of is the hydrophobic nature of the peptide seg- 
ments. In fact peptide insertions are usually discussed by 
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use of the hydrophobicity scale of the component amino 
acids [6, 7]. Other factors that control insertion such as 
membrane curvatures are also very intriguing [8, 9]. We 
should have it in mind that the insertion process depends 
sensibly on the primary structure of the peptide such as 
the lengths of the hydrophilic and hydrophobic sequences 
as well as on the properties of the matrix lipid membrane 
such as surface charge distributions and their imbalance 
between extracellular and cytoplasmic sides. 

Formations of various complexes in systems of water, 
surfactants, and amphiphilic block-copolymers are also 
of great interest from other points of view. The block- 
copolymers are known to stabilize bilayer structures of 
the L a phase of lipids [10]; it is also found that small 
addition of the block-copolymer drastically increase the 
emulsification capacity of the surfactants [11]. Of great 
biological and pharmaceutical interest is the finding that 
hydrophilic polymers with hydrophobic anchor greatly 
improve the stability the conventional liposome opening 
up new possibilities in drug delivery systems [12]. 

Despite large academic interest and potential ap- 
plications, complex cooperative actions of peptides 
(block-copolymers) and lipids that take place in local 
nanometer-sized space are often beyond the reach of ex- 
periments. Molecular simulation is now emerging as a 
very promising tool to study such complicated molecular 
processes. Molecular simulations of membrane proteins 
especially of their equilibrium structure and dynamics 
supporting biological functions have been investigated 
intensively [13], but there are only few reports on the 
molecular simulation studies of peptide insertion lead- 
ing finally to the folded structures [14-17]. Successful 



works were made mostly by modeling the membrane as 
a continuum or an aggregation of rigid cylindrical rods, 
in which detailed cooperative motions of the lipids and 
the peptide were not questioned [14-16]. 

We here consider the lipid molecule with a water- 
like hydrophilic head group and a polymethylene-like hy- 
drophilic tail, while the peptide is modeled as a block- 
copolymer made of hydrophobic and hydrophilic blocks 
of discrete hydrophobicity scales, though real peptides 
are made of amino acids having wide spectrum of hy- 
drophobicity. Once such an approach is found successful, 
however, the elaboration of the model into more realistic 
ones may be straightforward. 
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II. MODEL AND METHOD 

We will investigate the molecular mechanism of in- 
sertion of a model peptide into lipid bilayers by use of 
molecular dynamics (MD) simulation. We here adopt 
psudoatomic models for water, lipids, and peptide; char- 
acteristic length scale is of atomic dimension but the 
molecular structures are considerably simplified to ac- 
celerate computation. Lipid molecules that constitute 
the cell membrane are amphiphilic molecules that have 
both hydrophilic head groups and the hydrophobic tails. 
We consider a model lipid with a head group made of 
three hydrophilic beads (P-beads), and a tail made of 
seven hydrophobic beads (H-beads) (Fig. la). On the 
other hand the peptide is modeled as a block-copolymer 
made of hydrophilic blocks of the P-beads and hydropho- 
bic blocks of the H-beads (Fig. lb). The block-copolymer 
is assumed to be made of 200 beads and have a sequence 
P5(Hi5Pi5)eHi5; the length of the hydrophobic blocks of 
15 H-beads is so selected to give nearly the same length as 
the thickness of the bilayer in order to facilitate incorpo- 
ration into the membrane. The H-bead is here considered 
to be the usual united atom CH2 of mass 14, while the 
P-bead has also the same mass and interactions as the 
H-bead with only difference in the affinity to water. Each 
water molecule is represented by a single bead (W-bead) 
which has much larger interaction energy as described 
below. 

The lipid and peptide molecules are assumed to be 
fully flexible, where following bond-stretching potential 
is only imposed neglecting the energies of distortions in 
the bond angles and dihedral angles, 



Ubond = 7^( r M+l _ r o) 2 



(i) 



where r^+i is a length of the i-th bond, ro is the equi- 
librium bond length, and k^ is the spring constant. Fol- 
lowing non-bonded interactions are considered between 
two beads of the same chain more than two bonds apart, 
and between those of different chains. Between the beads 
of the similar hydrophobicity (P-P, H-H, W-W, and P- 
W), we assume the conventional Lennard- Jones interac- 
tions of the type, 



FIG. 1: (a) The model lipid is made of a head-group (three 
hydrophilic atoms P3: grey beads) and a tail (seven hy- 
drophobic atoms H7: black beads), (b) The model peptide 
is a copolymer made of hydrophilic (grey) and hydrophobic 
(black) blocks: P 5 (Hi 5 Pi5)6Hi 5 



TABLE I: Molecular parameters used in the simulation 



parameter 



in absolute units in reduced units 



&lipid — ^peptide 
O water 

Clipid — ^peptide 

twater 

771 



0.152 (nm) 
3.46 x 10 7 (J/mol/nm 2 ) 
0.38 (nm) 
0.33 (nm) 
500 (J/mol) 
2000 (J/mol) 
14.0 x 10" 3 (kg/mol) 



0.4 
9993 
1 

0.88 
1 

4.0 
1 



dW 



(2) 



, while between the hydrophilic and hydrophobic beads 
(H-P, and H-W) we impose soft core repulsion 



(3) 



Van der Waals attraction is here cut-off at a rather long 
distance as to assure the peptide's approach to the mem- 
brane. The water molecules are modeled as beads but 
with four times larger attractive interaction e in equation 
(2) considering approximately four time higher boiling 
point of water than that expected for hypothetical CH2 
liquid; the hydrophilic head of the lipid made of three 
beads is so constructed to have nearly the same total 
interaction with the water molecule. Present molecular 
model of lipid is largely based on a paper by Goetz and 
Lipowski [18, 19], but the molecular model is slightly fine 
grained to assure better correspondence to the usual hy- 
drocarbon chain which we have hitherto been studying. 
The potential parameters used in the present simulation 
is given in tabled I. 
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We made MD simulations under NTV ensemble with 
MD cell size 20<j in each direction and periodic bound- 
ary condition. We placed within the MD cell 305 lipids 
each made of 10 beads and a peptide molecule of 200 
beads together with surrounding 3675 water molecules. 
We solved Newton's equations of motion by the leapfrog 
method with temperature control at 300K by velocity 
scaling. We first confirmed that our present lipid model 
supports spontaneous formation of a bilayer membrane 
in water. Then we placed a model peptide molecule near 
the membrane and monitored the process of adsorption 
and insertion into the membrane. 



III. RESULTS AND DISCUSSIONS 
A. Spontaneous membrane formation 

The mesoscopic structure formation of lipids is a result 
of subtle balance of interactions among the hydrophilic 
heads, the hydrophobic tails, and water molecules. De- 
pending on the density of lipids, PH of the solvent, etc. 
various morphologies appear: micelles, vesicles, mem- 
branes, etc [1]. We first examined whether our molecu- 
lar model supports the spontaneous formation of bilayer 
membrane in water. It is readily noticed that, due to 
the applied periodic boundary condition, the membrane 
formation is possible only when the natural lateral width 
of the membrane fits the size of the MD cell, otherwise 
they will form spherical micelles, cylindrical micelles, etc 
[18]. By proper choice of the MD cell size (20a) and the 
number of lipids (305 chains), we could confirm a grad- 
ual formation of MD-cell-spanning bilayers out of random 
mixture of lipids and waters. In all the simulations that 
follow, we adopted this favorable lipid number and cell 
size that allow the membrane formation. 

Typical snapshot of the bilayers (Fig. 2a) obtained af- 
ter 1.0 ns of simulation at 300K shows that the observed 
membrane has large dynamical disorder that makes the 
quick response to external perturbation possible. The 
density profiles of the atoms are given in Fig. 2b, where 
both hydrophilic and hydrophobic atoms show rather 
wide distributions but still with a definite hydrophobic 
core. 



B. Peptide conformation before partitioning 

A free amphiphilic block-copolymer in water is con- 
sidered to have specific isotropic structure. Before dis- 
cussing insertion into the membrane, we give a glimpse at 
the conformation of our model peptide in water. Shown 
in Fig. 3 is a typical equilibrium conformation of the 
model peptide made of sequential hydrophilic and hy- 
drophobic blocks. The peptide is stretching six hy- 
drophilic block loops out into the water phase, while 
the hydrophobic blocks are segregated and shrunk to 
form spherical core. The amphiphilic peptide molecule is 




FIG. 2: (a) A snapshot of the lipid bilayers in water spon- 
taneously formed from a random mixture of water and lipid 
molecules at 300K. The z axis is taken perpendicular to the 
membrane, (b) The reduced number density (average num- 
ber of atoms per unit volume ), projected along the x-y plane, 
of the hydrophilic atoms (dash-dot) and hydrophobic atoms 
(solid) together with water (dash); the abscissa is also reduced 
by unit. 



thus imagined as a unimolecular micelle with hydropho- 
bic core segments surrounded by hydrophilic loops. This 
characteristic initial conformation of our model peptide 
in water has been profoundly distorted when it comes 
close to touch the bilayer membrane. The isotropic distri- 
bution of the hydrophilic loops around the hydrophobic 
core is disrupted, and the hydrophobic core gets bared of 
the hydrophilic overcoat resulting in strong interactions 
with the hydrophobic tails of the membrane. This also 
gives larger perturbation to the integrity and the pla- 
nar ity of the membrane, which finally leads to a success- 
ful translocation of the hydrophilic blocks of the peptide 
against the hydrophobic barrier in the membrane core as 
described below (Fig. 4). 
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FIG. 3: Typical chain conformation of the model peptide 
free from interactions with the membrane. The molecule is 
deemed as a unimolecular micelle stretching seven hydrophilic 
arms out into the water phase while the hydrophobic blocks 
are cramped to form a central core region. 



C. Peptide insertion into the membrane 

When the peptide molecule is set free to interact with 
the membrane, the interaction gives strong perturbation 
to the structures of both peptide and membrane near the 
point of contact. The hydrophobic blocks of the peptide 
strongly favor intimate contact with the lipid tails when 
they come closer than the interaction cut-off 4<r, while the 
attraction between the hydrophilic blocks of the peptide 
and the lipid heads is moderated by the intervening water 
molecules. 

The core of the unimolecular micelle peptide attracts 
lipid tails and tends to form even larger micelle result- 
ing in depletion of the lipid molecules around the con- 
tact point (Fig. 4c). Such large perturbation to the mem- 
brane often gives rise to the creation of a pore, though 
short-lived, in the membrane. Configuration of the chains 
around the pore is seen in Fig.4d, where lipid molecules 
tend to lie horizontally with their hydrophilic heads 
pointing outward into the water phase. We can notice 
that the side-surface of the pore is almost covered with 
the hydrophilic heads of the lipids, and this is considered 
to help the hydrophilic blocks of the peptide translocate 
across the membrane. The opening of the pore is only 
temporary and the membrane integrity is readily recov- 
ered in a very short time (Fig.4e). Finally the peptide 
acquires its stable position forming neat folded confor- 
mation within the membrane (Fig.4f). Figure 4 shows a 
whole trajectory of the peptide and the membrane dur- 
ing a successful insertion, while Fig. 5 shows the center-of- 
mass motion of the peptide. Around 0.2 ns, the peptide 
begins to contact the membrane and starts insertion, un- 
til it finally comes to stay in the middle of the membrane. 

Typical peptide conformation after 1.0 ns of simulation 
is given in Fig. 6a. The hydrophilic loops are translo- 




FIG. 4: Typical trajectories at 300K of the peptide success- 
fully translocated the bilayers: (a) the initial state at 0.1 ns, 
(b) start contacting at 0.2 ns, (c) uptake of lipid tails at 0.3 
ns, (d) pore formation at 0.7 ns, (e) recovering membrane 
integrity at 0.8 ns, (f) complete penetration at 1.7 ns. 



cated to the other side across the membrane and the 
hydrophobic blocks are nearly stretched perpendicular 
to the membrane, which is a typical conformation of the 
membrane peptide. The hydrophilic atoms of the peptide 
show nearly symmetric distribution (Fig. 6b), which in- 
dicates that about half of the hydrophilic blocks, two out 
of six in this case, are translocated across the membrane. 
We should here remark that the number of blocks that 
successfully penetrate the membrane is not always half 
the total hydrophilic blocks. Sometime only a small frac- 
tion of the blocks translocate the membrane resulting in 
quite asymmetric distribution. It is even the case when 
the peptide only adheres to the one side of the membrane 
and shows no sign of translocation. 

Lipid chains in the bilayers are known to show ac- 
tive diffusion within the layer, and with much less fre- 
quency between the layers; the later mode of chain dif- 
fusion is accompanied by the inversion of the chain head 
and ta il and it is called flip-flop motion. The fre- 
quency of the flip-flop motion is usually very low, but 
recent fluorescence experiment has revealed that it is 
greatly enhanced during the process of peptide insertion 
[20, 21]. We here monitored the rate of chain inversion 
(flip-flop); we counted thefrequency of the chain inver- 
sion that the head-to-tail vectors of the chains change 
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FIG. 5: Time evolution of the reduced z-coordinate of the 
center-of-mass of the peptide. Also shown are the positions 
of membrane surfaces (horizontal lines), which are represented 
by the centers of the hydrophilic groups of the lipids. 



FIG. 7: Total frequency of the flip-flop motions of the lipids 
within 2 ps interval; the data was averaged over ±50 ps in- 
terval. Marked increases are noticed during the contact and 
insertion of the peptide. 
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FIG. 6: (a) Typical peptide conformation that has succeeded 
in the translocation, and (b) distributions of the component 
hydrophilic (o) and hydrophobic (•) atoms along the z-axis. 



the sign of their z-component perpendicular to the mem- 
brane. Fig. 7 clearly shows that the flip-flop motion shows 
marked increases at the onset of peptide contact and in- 
sertion around 0.2 ns. Most of the chain inversion is cou- 
pled withthe formation of the pore, the surface of which 
is covered with lipids pointing their hydrophilic heads 
perpendicularly toward the water phase within the pore. 
The driving force for the insertion and translocation of 
peptides has in general many origins, and most of them 
require metabolic energy for specific machinery [2]. Even 
the short peptides which insert spontaneously can be 
driven by both energy and entropy. We here monitored 
the energy during the insertion process and decomposed 
it into characteristic interaction energies: interactions 
between (b) lipid tails and protein H-beads@^ a ^_i^, 
(c) lipid heads and protein P-beads Eh ea d-p, (d) lipid 
tails Etaii-taiU ( e ) protein H-beads E H -h, (f) protein P- 
beads and water Ep- water , as well as (a) total interaction 
energy E to tai (Fig. 8). Most conspicuous is the decrease 
in Etaii-H by the insertion of the peptide's hydropho- 
bic block into the hydrophobic membrane core. This en- 
ergy decrease is, however, counterbalanced both by the 
increase in E ta ii-taii (during the initial insertion up to 
0.5 ns) which indicates the disordering in the membrane 
core region, and by the increase in Eh-h (during the late 
stage of insertion after 0.5 ns) which is due to a slightly 
relaxed chain packing of the hydrophobic blocks of the 
peptide in the hydrophobic environment in the membrane 
core. The overall energy change is only slightly negative, 
and this suggests the possible importance of the entropic 
contribution to the insertionprocess. 

Taking all these discussions into account, the mecha- 
nism of spontaneous insertion of short peptide-like poly- 
mers can be described as follows. A free amphiphilic 
block-copolymer in water forms a unimolecular micelle 
with a collapsed hydrophilic core and extended hy- 
drophilic arms, though real membrane peptides are know 
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FIG. 8: Energy changes during the insertion process, and the decomposition into energies of different types of interactions; 
(a) total interaction energy , (b) the interaction between lipid tails and protein H-beads , (c) between lipid heads and protein 
P-beads , (d) between lipid tails , (e) between protein H-beads , and (f) between protein P-beads and water . The energy scale 
is the reduced unit, which is obtained by divided by . 



to be difficult to maintain such states in water due to 
their high hydrophobicity of the a-helix forming stems. 
When the polymer approaches the bilayers, flexible hy- 
drophilic arms give way to the hydrophobic core in con- 
tacting the membrane. Strong interactions between the 
core and the lipid tails perturb the membrane integrity 
resulting in a pore formation. The hydrophobic barrier to 
the peptide's translocation across the membrane is thus 
alleviated, and some of the hydrophilic blocks succeed in 
penetrating the membrane. Such molecular picture may 
be a little bit far from conventional ones obtained from 
other molecular simulations. But nice correspondence 
with recent experiments seems to support the present 
scenario. Furthermore, it is quite reasonable to consider 



that the amphiphilic lipid molecules of the membrane act 
as surfactants for the hydrophobic core of the peptide 
forming larger micelle. Thus the lipids and the peptide 
are required to perform highly cooperative actions during 
initial insertion, and through considerable reorganization 
of the lipids-peptide complex the peptide molecule com- 
pletes the insertion into the membrane. 
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